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(57) A method of estimating phase noise spectral 
density from a jitter versus time vector array obtained 
from a periodic signal having an average frequency in- 
cludes the an initial step of converting the jitter verus 
time vector array to a phase error versus time vector 
array using an estimate of the average frequency of the 
periodic signal. A time to frequency transform is applied 
to the phase error versus time vector array to generate 
a phase error magnitude versus frequency vector array, 
and a phase noise spectral density vector array ob- 
tained by normalizing the phase error magnitude versus 
time vector array to a one hertz bandwidth. The method 
of estimating the jitter in the periodic signal includes the 
steps of generating a vector array of estimated refer- 



ence crossing times of the periodic signal in the wave- 
form record using interpolation and calculating an esti- 
mated periodic signal frequency based on a selected 
number of reference crossings in the estimated refer- 
ence crossing times vector array and associated time 
positions of the waveform samples from the first select- 
ed reference crossing to the last selected reference 
crossing in the waveform record. A vector array of uni- 
formly spaced ideal reference crossing times is gener- 
ated based on the estimated periodic signal frequency, 
and a uniformly spaced jitter versus time vector array is 
generated by determining the difference between the 
ideal reference crossing times and the corresponding 
estimated reference crossing times of the periodic sig- 
nal. 
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Description 

CROSS-REFERENCE TO RELATED APPLICATION 

[0001] This application claims the benefit of the U.S. Provisional Application No. 60/1 77,750 filed January 20, 2000. 
BACKGROUND OF THE INVENTION 

[0002] The present invention related generally to the measurement of phase noise spectral density and jitter in a 
periodic signal and more particularly to the measurement of phase noise spectral density and jitter in a digitally sampled 
waveform record acquired in the time domain. 

[0003] Jitter is the time deviation of a periodic signal's "significant instants" from an ideal non-jittered periodic signal. 
The periodic signal may be an electrical or optical modulated or unmodulated carrier signal. Phase noise spectral 
density is the spectral density of unwanted phase modulation or jitter on the fundamental frequency of a periodic signal. 
The amount of jitter and phase noise in a periodic signal is an important measurement parameter for circuit and system 
designers. There are a number of measurement instruments and methods for measuring the jitter and phase noise. A 
spectrum analyzer measures the magnitude of a signal under test in the frequency domain. One drawback to using a 
spectrum analyzer to measure phase noise is that the instrument cannot distinguish between amplitude modulation 
noise and phase modulation noise. Additionally, conventional spectrum analyzers require a separate calibration pro- 
cedure to relate voltage or power to phase. It also has a small analysis frequency span around the carrier signal and 
the spectrum is not normalized to dB/carrier (dBc/Hz). 

[0004] Another system for measuring the phase error is a hardware implemented analog system having an accurate 
reference oscillator that outputs a periodic signal with the same frequency as the periodic signal being measured. The 
output of the oscillator is coupled to a multiplier circuit that receives the signal under test containing phase noise. The 
output of the multiplier is an signal having two components, the sum and difference of the reference oscillator output 
and the signal under test. The output of the multiplier is passed through a low pass filter to remove the high frequency 
components which leaves the phase noise of the signal under test. A drawback to this type of system is the requirement 
for the reference oscillator to have a stable and precise output. This would require a very precisely controlled phase 
locked loop (PLL). There would then be a need to correct the measured data for the PLL response for the offset 
frequencies less than the lock bandwidth. 

[0005] A further system employs a heterodyne/ counter method. A counter is triggered on a rising or falling edge of 
the periodic signal under test crossing a threshold value and incremented by an internal clock. On the next rising or 
falling edge crossing the threshold, the counter is stopped and the period of the internal lock times the count on the 
counter. On a subsequent rising or falling edge, the counter is again initiated and a count is acquired to the next rising 
or falling edge. One drawback to this type of system is that the results represent period disturbance. In addition, the 
system does not have a reference clock equal to the average frequency of the periodic signal under test. Further, the 
system does not produce uniformly space time intervals for the measured periods which is a necessity for accurate 
transform from the time domain to the frequency domain. Additionally, the system does not measure the time between 
each rising edge since each second edge is used to disable the counter. 

[0006] Digital oscilloscopes have also been used to measure the jitter of a signal under test in the time domain. An 
acquisition system acquires a waveform record of the signal under test. An average period of the signal under test is 
generated from the waveform record to produce an ideal clock. A reference crossing point is defined for the signal 
under test and the ideal clock. The magnitude difference between the reference crossings of the signal under test and 
reference crossings of the ideal clock are computed to produce magnitude phase error values relative to the signal 
under test. The magnitude phase error values need to be resampled at an interpolated sample rate to generate an 
equally spaced sample interval that is applied to a time to frequency transform to generate a magnitude versus fre- 
quency output of the magnitude phase error values. A drawback to the above described method is the need to resample 
the error data at an interpolated sample rate to generate uniformly space phase error values. Further, the signal under 
test may contain amplitude noise and additive noise from the measurement instrument in addition to phase noise. Such 
instrument generated noise cannot be separated from the phase noise and will be transformed as part of the phase 
noise spectral density. 

[0007] What is needed is an improved method for measuring the phase noise spectral density of a periodic signal 
that does not require the interpolated resampling of the error values to produce uniformly spaced error values. Addi- 
tionally, the improved method should provide a new method of determining the jitter in the periodic signal that is not 
dependent on measured magnitude values. 
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SUMMARY OF THE INVENTION 

[0008] Accordingly, the present invention is a method of estimating phase noise spectral density and jitter in a periodic 
signal acquired in a digitally sampled waveform record in the time domain. The present invention estimates the phase 
noise spectra! density and jitter by relying on timing information of events in the digitally sampled waveform record. 
One aspect of the invention is estimating the jitter in the periodic signal that includes the step of generating a vector 
array of estimated reference crossing times of the periodic signal in the waveform record using interpolation. An esti- 
mate of the periodic signal frequency is calculated based on a selected number of reference crossings in the estimated 
reference crossing times vector array and associated time positions of the waveform samples from the first selected 
reference crossing to the last selected reference crossing in the waveform record. A vector array of uniformly spaced 
ideal reference crossing times is generated based on the estimated periodic signal frequency, and a uniformly spaced 
jitter versus time vector array is generated by determining the difference between the ideal reference crossing times 
and the corresponding estimated reference crossing times of the periodic signal. 

[0009] A further aspect of the invention is estimating phase noise spectral density from the jitter versus time vector 
array that includes the step of converting the jitter verus time vector array to a phase error versus time vector array 
using an estimate of the average frequency of the periodic signal in radians per second. A time to frequency transform 
is applied to the phase error versus time vector array to generate a phase error magnitude versus frequency vector 
array, and a phase noise spectral density vector array is obtained by normalizing the phase error magnitude versus 
time vector array to a one hertz bandwidth. A frequency band may be defined in the phase noise spectral density vector 
array, and the phase noise values in the phase noise spectral density vector array within the defined frequency band 
are integrated to obtain a jitter RMS value within the defined frequency band. 

[0010] The generation of the estimated reference crossing times vector array step further includes the alternative 
steps of either estimating rising edge reference crossing times of the periodic signal in the waveform record or falling 
edge reference crossing times of the periodic signal in the waveform record. A further alternative step includes con- 
currently generating a vector array of estimated rising reference crossing times of the periodic signal in the waveform 
record using interpolation and a vector array of estimated falling reference crossing times of the periodic signal in the 
waveform record using interpolation. From the rising and falling reference crossing times vector arrays a vector array 
of estimated pulse width deviation times is generated by comparing the estimated rising reference crossing times vector 
array to the estimated falling reference crossing times vector array. 

[001 1 ] The interpolation to generated the estimated reference crossing times may be implemented by linearly inter- 
polating between the first data sample above the reference crossing and the first data sample below the reference 
crossing. Alternately, the interpolating step may be implemented by generating a higher order interpolator using multiple 
digital data samples above and below the reference crossing where the higher order interpolator step includes the step 
of generating a window, such as a sin(x)/x function. 

[001 2] The calculating of the estimated periodic signal frequency includes the step of interpolating an average slope 
from the number of reference crossing times and the sum of the reference crossing times where the interpolation step 
includes calculating a best fit linear curve to the respective number of reference crossing times and the corresponding 
reference crossing times, and estimating the periodic signal frequency from the slope of the fitted line. 
[001 3] The objects, advantages and novel features of the present invention are apparent from the following detailed 
description when read in conjunction with appended claims and attached drawings. 

BRIEF DESCRIPTION OF THE DRAWING FIGURES 

[0014] Figure 1 is a representative block diagram of a time domain measurement system implementing the method 
of estimating the phase noise spectral density and jitter in a periodic signal according to the present invention. 
[0015] Figure 2 is a representative waveform trace of a periodic signal having phase errors. 

[0016] Figure 3 is an enlarged view of a rising edge of the periodic signal crossing a reference level illustrating the 
interpolation of the crossing point time of the rising edge. 

[0017] Figure 4 is a graphical representation for estimating the frequency of the periodic signal. 
[0018] Figure 5 is a graphical representation of the generated uniformly spaced ideal reference crossings from the 
estimated periodic signal frequency and the jitter between the ideal reference crossings and the estimated reference 
crossings of the periodic signal. 

[0019] Figures 6 is a graphical representation of the jitter versus time of the uniformly spaced ideal reference cross- 
ings. 

[0020] Figures 7A and &B are a flow chart showing the steps in estimating the phase noise spectral density and jitter 
in a periodic signal according to the present invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0021] In the following detailed description numerous specific details are set forth in order to provide a thorough 
understanding of the present invention. However, it will be understood by those skilled in the art that the present 
invention may be practiced without these specific details. In other instances well known methods, procedures, com- 
ponents, and circuits have not been described in detail so as not to obscure the present invention. 
[0022] Part of the description will be presented in terms of operations performed by a measurement instrument, using 
terms such as data, values, signal samples, numbers and the like, consistent with the manner commonly employed 
by those skilled in the art to convey the substance of their work to others skilled in the art. As is well understood by 
those skilled in the art, these quantities take the form of electrical, magnetic, or optical signals capable of being stored, 
transferred, combined and otherwise manipulated through mechanical and electrical components of the measurement 
instrument; and the term measurement instrument includes general purpose as well as special purpose data processing 
machines, systems, and the like, that are stand alone, adjunct or embedded. 

[0023] Additionally, various operations will be described as multiple discrete steps in turn in a manner that is most 
helpful in understanding the present invention, however, the order of description should not be construed as to imply 
that these operations are necessarily order dependent, in particular, the order of their presentation. 
[0024] Referring to Figure 1, there is shown a representative block diagram of a digital sampling system 10 that 
acquires a waveform record of an input periodic signal under test 12. The measurement system may be a digital 
oscilloscope having an acquisition subsystem, representatively shown at 14, and an instrument control subsystem 16. 
Alternately, the acquisition subsystem 14 may be a separate digitizing unit with the instrument control subsystem 16 
being a host computer that receives acquired waveform data from the digitizer. The input periodic signal 12 is coupled 
through a variable attenuator 18 and a preamplifier 20. In high digitizing rate sampling oscilloscopes, such as the 
TDS71 04 Digital Oscilloscope, manufactured and sold by Tektronix, Inc. Beaverton, Oregon and assignee of the instant 
invention, each input channel has digitizing pipes, as representatively shown as pipes 22 and 24. Any number of pipes 
may be included for each oscilloscope input channel. Each pipe has a track-and-hold (T/H) circuit 26, an analog-to- 
digital (A/D) converter 28 and a memory 30. A time base 32 provides timing signals to the acquisition subsystem 14 
for latching an analog value of the input signal in the T/H circuits 26, clocking the A/D converters 28 to digitize the 
analog value on the T/H circuit and storing the digitized values in memories 30. The pipes 22, 24 have additional 
circuitry (not shown) that offsets or delays the time base signals for each pipe 22, 24 to produce an acquisition rate 
"X" times the timing signal rate where "X" is the number of pipes in the acquisition subsystem 14. 
[0025] The time base 32 is coupled to a control and address bus 34 for receiving commands from a DSP controller 
36. The DSP controller may be an application specific integrated circuit (ASIC) or a microprocessor, such as a Intel 
CELERON™ microprocessor, manufactured and sold by Intel, Corp., Santa Clara, CA. Also coupled to the control and 
address bus is the variable attenuator 18 and the circuitry within the digitizing pipes 22, 24. The acquisition subsystem 
14 may also include acquisition controls 38 and memory 40 when the subsystem 14 is a separate digitizing unit. The 
acquisition controls 38 may include buttons, rotatable knobs and the like and/or control entry devices, such as a key- 
board and/or mouse, for setting the time base and acquisition parameters for the digitizing unit. The memory 40 contains 
stored program instructions that are accessed by the DSP controller 36 for controlling the operation of the acquisition 
subsystem 14. When the acquisition subsystem 14 is part of a digital oscilloscope, the acquisition controls 38 and 
memory 40 are incorporated into the front panel and system memory of the oscilloscope. An external reference line 
42 is coupled to the time base 32 that allows the time base to be locked to an external reference source. 
[0026] The acquisition subsystem 1 4 is coupled to the instrument control subsystem 1 6 via an interface bus 44. The 
interface bus 44 provided bi-directional communications for transferring data and control signal between the subsys- 
tems. The interface bus 44 may be a commercially available bus adapter integrated circuit that controls the transfer of 
data between two bus systems. The interface bus may also be a bidirectional serial bus, such as an internal l 2 C bus 
for an oscilloscope or an external GPIB bus for a digitizing unit and a computer. The instrument control subsystem 1 6 
includes a system or host controller 46, coupled via system bus 48 to system memory 50, a display device 52, front 
panel controls 54, and a mass storage unit 56. The system memory 50 includes both RAM, ROM and cache memory 
with the RAM memory storing volatile data, such as the data values representative of the input signal passed by the 
acquisition subsystem 14. The controller 46 is preferably a microprocessor, such as PENTIUM® microprocessor, man- 
ufactured and sold by Intel, Corp., Santa Clara, CA. The display device 52 may be a liquid crystal display, cathode ray 
tube or the like, and front panel controls 54 may include buttons, rotatable knobs and the like and/or control entry 
devices, such as a keyboard and/or mouse. The mass storage unit or units 56 may be a hard disk drive, a CD ROM 
drive, a tape drive, a floppy drive or the like, that reads from and/or writes to appropriate mass storage media. The 
digital sampling phase error measurement system 10 in the preferred embodiment of the invention is a PC based 
system controlled under WINDOWS® 98 operating system, manufactured and sold by Microsoft, Corp., Redmond, WA. 
[0027] The present invention has the distinct advantage of estimating the phase noise spectral density and jitter by 
relying on timing information of events in the digitally sampled waveform record. Additionally, the present invention is 
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a computationally efficient of producing uniformly spaced jitter and phase noise spectral density data using the time 
data of the digitally acquired waveform record. 

[0028] The below description of the method of estimating the phase noise spectral density and jitter in a digitally 
acquired waveform record of a periodic signal will be described with reference to Figures 2 through 7. Referring to 

5 representative waveform trace of Figures 2, the is shown a periodic signal 70 having rising and falling edges 72 and 
74 representing a dipped sinusoidal signal The clipped sinusoidal signal is for illustrative purposed only and any type 
of periodic signal may be used in the below described method. Also for illustration purposes, the periodic signal is 
shown with a DC level. The periodic signal is referenced to a reference level, labeled 76, which may or may set to zero 
volts providing the periodic waveform is not at a DC level. The periodic signal 70 is sampled at a uniform sample rate 

10 by the acquisition subsystem 14 as represented by the "x"s on the waveform trace. The vertical dashed fines 77 rep- 
resent the reference crossings of the periodic signal 70 without jitter. In the preferred embodiment, the acquisition 
subsystem 14 can sample the periodic signal at up to a 20 Gsamples/sec rate. That would allow a 100 MHZ periodic 
signal to be sampled 200 times per period. Usually, the reference level 76 will fail between two successive waveform 
samples on the rising or falling edges of the periodic signal requiring interpolation to estimate the crossing point of the 

15 edge relative to the reference level. 

[0029] Figure 3 shows an enlarged view of a periodic signal edge 78 crossing the reference level 76. As can be seen 
in the figure, there is no waveform sample at the reference level. The simplest form of interpolation is linear interpolation 
of at least a first data sample 80 above the reference crossing point and a first data sample 82 below the reference 
crossing point. The time difference between the data sample above the reference and the data sample below the 

20 reference level is calculated. A weighting factor "A" derived from the slope of the line between the data samples above 
and below the reference level is applied to the calculated time difference and the resulting time value is the estimated 
time of the reference level crossing of the periodic signal edge Higher order interpolation may be used for determining 
the time of the reference crossing point, such as a windowed sin(x)/x function using multiple samples above and below 
. the reference crossing point. 

25 [0030] The interpolated reference crossing times of the periodic signal edges {t 1 f eai, ....V 631 ) ar © stored in a vector 
array called tjret_real as represent in box 90 in Figure 7A. 

LreLreat = [t, rea ', t 2 ma \ t^*', ... t n reat ) (1) 

30 

The Lref_real vector array may be composed of either the reference crossing times of the rising edges of the period 
signal or the falling edges of the periodic signal. It is also possible to estimate the reference crossing levels for both 
the rising and falling edges of the periodic signal and concurrently generate estimated reference crossing times vector 
arrays for the rising edges t_rise_ref_real and the falling edges t_fall_ref_real. 



35 



45 



50 



(2) 



UalLreLreal = {<-. fi*,. G*. - O P> 

A pulsewidth deviation vector array t_pw_dev can be generate by comparing the rising edge reference crossing vector 
array t_rise_ref_reaf in relation to the falling edge reference crossing vector array t_fall_ref_reat. 

t_pw_dev = {% w - de \ f™~ de \ f™- dev , ... f n w - dev } {4 ) 

[0031] The method of estimating the periodic signal frequency from the t_ref_rea\ vector array is best understood 
with reference to Figure 4. The frequency of the periodic signal is estimated using the reference crossings in the 
Lref_real reference crossing vector array. If the t_ref_real reference crossing vector array is large indicating a large 
waveform record, a smaller number of reference crossing may be used so long as a sufficient number of reference 
crossings are selected to provide a statistically significant representation of the periodic signal period. The horizontal 
axis is scaled in reference crossing numbers and the vertical scaled in time. The reference level crossing times are 
plotted to the respective crossing number using a linear scale on both axes. As shown in the figure, the plot of the 
reference level crossing times reveals an approximate linear curve that represents the frequency of the periodic signal. 
A best fit linear curve is generated for the Cref_real vector array and the slope of the fitted line is the estimate of the 
periodic signal frequency. The linear curve may be generated by taking the number of periods in the selected time 



5 



EP1 118 866 A1 



interval and dividing by the time difference between the reference level crossing at t n and t 1 as represented by box 92 
in Figure 7A and the following equation. 



f _ 1 _ (nl) 



c ,jeal- seal.. ( 5 ) 



[0032] The calculated frequency of the periodic signal f c is used to generate a vector array of uniformly spaced ideal 
reference crossing points t_ref_ideal of the periodic signal. The initial reference crossing point t^ 3 ' of the periodic 
signal is used as the initial ideal crossing point. The vector array t_refjdeal contains the ideal time series that would 
correspond to phase error free uniformly spaced reference crossing as represented box 94 in Figure 7 A and the below 
vector sequence. 

LrefJdeal = {/f*', /f ea ' + + * ,f ea ' + 3 ... ,f + <JL_1>, (6) 

'c f c f c 

[0033] Figures 5 and 6 graphically show the generation of a uniformly spaced jitter versus time jitter_time vector 
array. Both the horizontal and vertical axis are scaled in time with the horizontal axis divided into uniformly space ideal 
reference crossing point from the t_ref_ideal vector array. The reference crossing times from the t_ref_real reference 
crossing vector array are plotted on the vertical axis. The diagonal line represents the frequency of the periodic signal. 
tjdeai \ s obtained as the time position or coordinate of the intersection of the estimated ideal frequency with the horizontal 
axis in Figure 5. The difference between the intersections of the P deal times and the corresponding V eai times and the 
reference frequency line is the jitter. The jitter^time vector array is uniformly spaced on the i deai times. The jitter_time 
vector array in the generated by subtracting on an element by element basis the periodic signal real reference crossing 
points from the corresponding ideal reference crossing points are represented by box 96 in Figure 7A and below vector 
sequence. 



(7) 



[0034] The jitter_time vector array is converted to a phase error versus time vector array phase_error__time by mul- 
tiplying the elements of the jitter_time vector array by the periodic signal frequency in radians per second as represented 
by box 98 in Figure 7 A and the below vector sequence. 

phase_error_time = { 9l ( t, idea ), . . . <p n ( t n ideat )} ( 8 ) 

where <p, = cp (f? eal ), i, n . The significance of the phase_error_time vector array is that it was uniformly spaced in time 
and generated from jitter_time vector array that is scaled in time and not in voltage magnitude. Having an uniformly 
spaced phase_error_time vector array allows the generation of an accurate frequency domain spectral density vector 
array spectrat_phase_noise from the phase_error_time vector array. A time to frequency transform, such as a Fast 
Fourier Transform function, is applied to the phase_error_time vector array to generate a phase error magnitude versus 
frequency vector array phase_error_freq as represented by box 1 00 in Figure 7A and the below vector sequence: 

phase_error_freq = {<p(F,), . . . <p(F„)} (9) 

2 

where F, = (i - 1) x F c /2. There are a number of time to frequency transforms that can convert a time domain sequence 
to the frequency domain sequence and the method of the present invention is not limited to any one particular transform. 
[0035] The phase_error_freq vector array is normalized to the to a 1 Hz bandwidth to generate a spectral density 
phase noise vector array spectraLphase^noise, as represented by box 102 in Figure 7A and the below vector se- 
quences 



spectraLphase^oise^ {<pF,- x ~ 1 ^ } (io) 
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where / denotes frequency components in vector sequence (9) that are not deterministic (i.e. a pure tone). 

spectral_phase_noise j = (q> F { ) i (11) 

where / denotes deterministic frequency components in vector sequence (9) (i.e. discrete frequency components). 
Those skilled in the art might employ different methods for identifying discrete frequency components in vector se- 
quence (9). One such method would sub-sample the phase_error_time vector array, re-compute the 
spectral_phase_noise vector array using only the formula in vector sequence (10) and identifying frequency compo- 
nents in the spectral_phase_noi$e vector array which remain constant as discrete components. 
[0036] Jitter RMS may be obtained within a frequency band of phase noise spectral density vector array by using a 
filter function to define a frequency band within the spectral _phase_noise vector array as represented by box t04 in 
Figure 7B and performing a numerical integration on the filtered phase noise values over the defined frequency band 
as represented by box 106 in Figure 7B. Numerical integration can also be performed on the all the elements in the 
spectral_phase_noise vector array to obtain the total jitter RMS of the periodic signal. 

[0037] A method has been described for estimating phase noise spectral density and jitter in a digitally sampled 
periodic signal. The method includes the steps of generating a vector array of the estimated reference crossing times 
of the periodic signal in the waveform record and calculating an estimated periodic signal frequency based on a esti- 
mated reference crossing time of the periodic signal. The estimated periodic signal frequency is used to generate a 
vector array of uniformly spaced ideal crossing times. A uniformly spaced vector array of jitter versus time is generated 
by determining the difference between the ideal crossing times and the corresponding estimated reference crossing 
times of the periodic signal. The jitter versus time vector array is converted to a phase error versus time vector array 
by multiplying the element of the vector array by the frequency of the periodic signal in radians per second. A phase 
error magnitude versus frequency vector array is generated by applying a time to frequency transform function to the 
phase error versus time vector array. The phase error magnitude versus frequency vector array is normalized to a one 
hertz bandwidth to produce a phase nosie spectral density vector array. Numerical integration may be performed on 
ail or a portion of the phase noise spectral density vector array to generate RMS jitter values for the periodic signal or 
a frequency band within the phase noise spectral density vector array. 

[0038] It will be obvious to those having skill in the art that many changes may be made to the details of the above- 
described embodiments of this invention without departing from the underlying principles thereof. For example, the 
digital sampling system may be synchronously sample the periodic signal under test by locking the time base to an 
external reference source having a frequency equal to the periodic signal under test. A vector array of voltage versus 
estimated reference crossing times may then be generated the waveform record and a jitter versus time vector array 
generated by dividing the voltage versus time array by the slope of the periodic signal in the vicinity of the crossings. 
The scope of the present invention should, therefore, be determined only by the following claims. 



Claims 

1 . A method of estimating phase noise spectral density from a jitter versus time vector array obtained from a periodic 
signal having an average frequency comprising the steps of: 

a) converting the jitter versus time vector array to a phase error versus time vector array using an estimate of 
the average frequency of the periodic signal; 

b) applying a time to frequency transform to the phase error versus time vector array to generate a phase 
error magnitude versus frequency vector array; and 

c) obtaining a phase noise spectral density vector array by normalizing the random phase error magnitude 
component of the phase error magnitude versus frequency vector array to a one hertz bandwidth. 

2. The method of estimating phase noise spectral density as recited in claim 1 further comprising the steps of: 

a) using a filter function to define a frequency band in the phase noise spectral density vector array; and 

b) integrating the filtered phase noise values in the phase noise spectral density vector array over the defined 
frequency band to obtain a jitter RMS value within the define frequency band. 

3. The method of estimating phase noise spectral density as recited in claim 1 wherein the converting step further 
comprises the step of multiplying the jitter versus time vector array by average frequency of the periodic signal in 
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radians per second. 

4. A method of estimating phase noise spectral density in a periodic signal acquired in a digitally sampled waveform 
record comprising the steps of: 

a) generating a vector array of estimated reference crossing times of the periodic signal in the waveform record 
using interpolation; 

b) calculating an estimated periodic signal frequency based on a selected number of reference crossings in 
the estimated reference crossing times vector array and associated time positions of the waveform samples 
from the first selected reference crossing to the last selected reference crossing in the waveform record; 

c) generating a vector array of uniformly spaced ideal reference crossing times based on the estimated periodic 
signal frequency; 

d) generating a uniformly spaced jitter versus time vector array by determining the difference between the 
ideal reference crossing times and the corresponding estimated reference crossing times of the periodic signal; 

e) converting the jitter verus time vector array to a phase error versus time vector array using an estimate of 
the average frequency of the periodic signal; 

f) applying a time to frequency transform to the phase error versus time vector array to generate a phase error 
magnitude versus frequency vector array; and 

g) obtaining a phase noise spectral density vector array by normalizing the random phase error magnitude 
component of the phase error magnitude versus frequency vector array to a one hertz bandwidth. 

5. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 wherein the gen- 
erating of the estimated reference crossing times vector array step further comprises the step of estimating rising 
edge reference crossing times of the periodic signal in the waveform record. 

6. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 wherein the gen- 
erating of the estimated reference crossing times vector array step further comprises the step of estimating falling 
edge reference crossing times of the periodic signal in the waveform record. 

7. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 wherein the gen- 
erating of the estimated reference crossing times vector array step further comprises the steps of: 

a) concurrently generating a vector array of estimated rising reference crossing times of the periodic signal in 
the waveform record using interpolation and a vector array of estimated falling reference crossing times of the 
periodic signal in the waveform record using interpolation; and 

b) generating a vector array of estimated pulse width deviation times by comparing the estimated rising ref- 
erence crossing times vector array to the estimated falling reference crossing times vector array. 

8. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 wherein the gen- 
erating of the estimated reference crossing times vector array step further comprises the step of interpolating 
between at least a first data sample above the reference crossing and at least a first data sample below the ref- 
erence crossing. 

9. The method of estimating phase noise spectral density in a periodic signal as recited in claim 8 wherein interpolating 
step further comprises the step of linearly interpolating between the first data sample above the reference crossing 
and the first data sample below the reference crossing. 

10. The method of estimating phase noise spectral density in a periodic signal as recited in claim 8 wherein the inter- 
polating step further comprises the step of generating a higher order interpolator using multiple digital data samples 
above and below the reference crossing. 

11. The method of estimating phase noise spectral density in a periodic signal as recited in claim 10 wherein the 
generating a higher order interpolator step further comprises the step of generating a windowed sin(x)/x function. 

12. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 wherein the cal- 
culating the estimated periodic signal frequency further comprising the step of interpolating an average slope from 
the number of reference crossing times and the sum of the reference crossing times. 
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13. The method of estimating phase noise spectral density in a periodic signal as recited in claim 12 wherein the 
interpolating step further comprises the steps of: 

a) calculating a best fit linear curve to the respective number of reference crossing times and the corresponding 
reference crossing times; and 

b) estimating the periodic signal frequency from the slope of the fitted line. 

14. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 further comprising 
the steps of: 

a) using a filter function to define a frequency band in the phase noise spectral density vector array; and 

b) integrating the filtered phase noise values in the phase noise spectral density vector array over the defined 
frequency band to obtain a jitter RMS value within the define frequency band. 

15. The method of estimating phase noise spectral density in a periodic signal as recited in claim 4 wherein the con- 
verting step further comprises the step of multiply the jitter versus time vector array by average frequency of the 
periodic signal in radians per second. 

16. A method of estimating jitter in a periodic signal acquired in a digitally sampled waveform record comprising the 
steps of: 

a) generating a vector array of estimated reference crossing times of the periodic signal in the waveform record 
using interpolation; 

b) calculating an estimated periodic signal frequency based on a selected number of reference crossings in 
the estimated reference crossing times vector array and associated time positions of the waveform samples 
from the first selected reference crossing to the last selected reference crossing in the waveform record; 

c) generating a vector array of uniformly spaced ideal reference crossing times based on the estimated periodic 
signal frequency; and 

d) generating a uniformly spaced jitter versus time vector array by determining the difference between the 
ideal reference crossing times and the corresponding estimated reference crossing times of the periodic signal; 

17. The method of estimating jitter in a periodic signal as recited in claim 1 6 wherein the generating of the estimated 
reference crossing times vector array step further comprises the step of estimating rising edge reference crossing 
times of the periodic signal in the waveform record. 

18. The method of estimating jitter in a periodic signal as recited in claim 16 wherein the generating of the estimated 
reference crossing times vector array step further comprises the step of estimating falling edge reference crossing 
times of the periodic signal in the waveform record. 

19. The method of estimating jitter in a periodic signal as recited in claim 16 wherein the generating of the estimated 
reference crossing times vector array step further comprises the steps of: 

a) concurrently generating a vector array of estimated rising reference crossing times of the periodic signal in 
the waveform record using interpolation and a vector array of estimated falling reference crossing times of the 
periodic signal in the waveform record using interpolation; and 

b) generating a vector array of estimated pulse width deviation times by comparing the estimated rising ref- 
erence crossing times vector array to the estimated falling reference crossing times vector array. 

20. The method of estimating jitter in a periodic signal as recited in claim 16 wherein the generating of the estimated 
reference crossing times vector array step further comprises the step of interpolating between at least a first data 
sample above the reference crossing and at least a first data sample below the reference crossing. 

21 . The method of estimating jitter in a periodic signal as recited in claim 20 wherein interpolating step further comprises 
the step of linearly interpolating between the first data sample above the reference crossing and the first data 
sample below the reference crossing. 

22. The method of estimating jitter in a periodic signal as recited in claim 20 wherein the interpolating step further 
comprises the step of generating a higher order interpolator using multiple digital data samples above and below 
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the reference crossing. 

23. The method of estimating jitter in a periodic signal as recited in claim 22 wherein the generating a higher order 
interpolator step further comprises the step of generating a windowed sin(x)/x function. 

24. The method of estimating jitter in a periodic signal as recited in claim 20 wherein the calculating the estimated 
periodic signal frequency further comprising the step of interpolating an average slope from the number of reference 
crossing times and the sum of the reference crossing times. 

25. The method of estimating jitter in a periodic signal as recited in claim 24 wherein the interpolating step further 
comprises the steps of: 

a) calculating a best fit linear curve to the respective number of reference crossing times and the corresponding 
reference crossing times; and 

b) estimating the periodic signal frequency from the slope of the fitted line. 
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